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Abstract. In at least two black hole binary systems, the Rossi X-Ray Timing Explorer has detected high frequency quasi- 
periodic oscillations (HFQPOs) with a 2:3 frequency commensurability. We propose a simple hot spot model to explain the 
positions, amplitudes, and widths of the HFQPO peaks. Using the exact geodesic equations for the Kerr metric, we calculate 
the trajectories of massive test particles, which are treated as isotropic, monochromatic emitters in their rest frames. By varying 
the hot spot parameters, we are able to explain the different features observed in "Type A" and "Type B" QPOs from XTE 
J1550-564. In the context of this model, the observed power spectra allow us to infer values for the black hole mass and 
angular momentum, and also constrain the parameters of the model. 



INTRODUCTION 

One of the most exciting results from the Rossi X-Ray 
Timing Explorer was the discovery of high frequency 
quasi-periodic oscillations from neutron star and black 
hole binaries |1, 2, 3]. For black hole systems, these 
HFQPOs are observed repeatedly at more or less con- 
stant frequencies, and in a few cases with integer ratios 
0]. These measurements give the exciting prospect of 
determining a black hole's mass and spin, as well as test- 
ing general relativity in the strong-field regime. To quan- 
titatively interpret these observations in the context of 
the space-time behavior near the black hole inner-most 
stable circular orbit (ISCO), we have developed a ray- 
tracing code to model the X-ray light curve from a col- 
lection of "hot spots," small regions of excess emission 
moving on geodesic orbits. The methods and basic re- 
sults of this code are described in detail in ||5|. 

This hot spot model is motivated by the similarity be- 
tween the QPO frequencies and the black hole coordinate 
frequencies near the ISCO |6, 7] as well as the suggestion 
of a resonance leading to integer commensurabilities be- 
tween these coordinate frequencies |8, 9]. Stella and Vi- 
etri 1 7] investigated primarily the QPO frequency pairs 
found in LMXBs with a neutron star accretor, but their 
basic methods can be applied to black hole systems as 
well. Perhaps the most powerful feature of this hot spot 
model is the facility with which it can be developed and 
extended to more general accretion disk geometries. For 
example, our ray-tracing code could be used in conjunc- 
tion with a 3D MHD calculation of the accretion disk 
to simulate the time-dependent X-ray flux and spectrum 
from such a disk. 



RAY TRACING IN THE KERR METRIC 

While the details of the ray tracing algorithm are pre- 
sented in full in jsl, we give here a short summary of 
our methods. We begin by dividing the image plane into 
regularly spaced "pixels" of equal solid angle in the ob- 
server's frame, each corresponding to a single ray. Fol- 
lowing the sample rays backward in time, we calculate 
the original position and direction that a photon emit- 
ted from the disk would require in order to arrive at 
the appropriate position in the detector The gravitational 
lensing and magnification of emission from the plane 
of the accretion disk is performed automatically by the 
geodesic integration of these evenly spaced photon tra- 
jectories, so that high magnification occurs in regions 
where nearby points in the disk are projected to points 
with large separation in the image plane. To model the 
time-varying emission from the disk, each photon path 
is marked with the time delay along the path from the 
observer to the emission point in the disk. 

To integrate the geodesic trajectories of photons or 
massive particles, we use a Hamiltonian formalism that 
takes advantage of certain conserved quantities in the 
dynamics: the energy, angular momentum, and mass of 
each particle. By eliminating the energy {E = —pt) from 
the Hamiltonian, we can use the conjugate coordinate t as 
the integration variable, thus reducing the problem from 
eight phase space dimensions to six. Another integral of 
the motion. Carter's constant Q ifioll . is not explicitly 
conserved in our equations, but rather is used as an 
independent check for the accuracy of the numerics. 
The resulting equations of motion do not contain any 
sign ambiguities from turning points in the orbits, as are 
introduced by many classical treatments of the geodesic 
equations in the Kerr metric. 



The photon trajectories are integrated backward in 
time from the image plane oriented at some incHnation 
angle / with respect to the axis of rotation for the black 
hole, where / = 0° corresponds to a face-on view of the 
disk and / = 90° is an edge-on view. The accretion disk 
is confined to a finite region of latitude and is assumed 
to be oriented normal to the rotation axis. The photons 
terminate either at the event horizon or pass through 
the surfaces of colatitude 9 — const. As trajectories pass 
through the disk, the photon's position and momentum 
(x^ ,Plj,) are recorded for each plane intersection in order 
to later reconstruct an image of the disk. The results of 
this paper are based primarily on flat disks with simple 
hot spot perturbations, yet with the computational meth- 
ods described above, arbitrary disk geometries and emis- 
sivity/opacity models can be simulated as well. 



HOT SPOT EMISSION 

Given the map from the accretion disk to the image 
plane, with each photon bundle labeled with a distinct 
4-momentum and time delay, we can reconstruct time- 
varying images of the disk based on time-dependent 
emission models. The simplest model we consider is 
a single region of isotropic, monochromatic emission 
following a geodesic trajectory: the "hot spot" model 
13, V] . We treat the hot spot as a small region of the 
disk with additional emissivity chosen to have a Gaussian 
distribution in local Cartesian space. We typically take 
Rspot = 0.25 — 0.5M, but find the normalized light curves 
and QPO power spectra to be independent of spot size 
and shape. 

By definition the hot spot will have a higher tempera- 
ture or density and thus greater emissivity than the back- 
ground disk, adding a small modulation to the total flux. 
RXTE observations find the HFQPO X-ray modulations 
to have typical amplitudes of 1-5% of the mean flux dur- 
ing the outburst |4]. Assuming a Shakura-Sunyaev type 
disk 1 11] with steady-state emissivity g{r) °^ r^^ and a 
similar scaling for the hot spot emission, we find that 
small hot spots with overbrightness of sa 100% (consis- 
tent with 3D MHD calculations LlZl ) are easily capable 
of creating X-ray modulations on the order of 1% rms. 

As the disk inclination increases, the light curve goes 
from nearly sinusoidal to being sharply peaked by special 
relativistic beaming. Thus the shape of a QPO light curve 
may be used to determine the disk inclination. Since cur- 
rent observational cannot resolve the X-ray signal over 
individual periods as short as 3-5 ms, instead the Fourier 
power spectrum can be used to identify the harmonic 
features of the light curve over many orbits. Disks with 
higher inclinations will give more power in the higher 
harmonic frequencies, due to the "lighthouse" effect, as 
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FIGURE 1. Fourier amplitude a„(rms) in higher harmonic 
frequencies v„ = nVij, as a function of orbital inclination to 
the observer, normalized as in Eqn. |2| The hot spot has size 
^spot = 0.5M, an overbrightness factor of 100%, and is in a 
circular orbit at i?isco around a Schwarzschild black hole. 

the hot spot shoots a high-power beam of photons to- 
ward the observer once per orbit, approximating a peri- 
odic delta-function in time. 

Figure [2 shows the quantitative dependence of har- 
monic power on disk inclination for a hot spot orbiting 
at the ISCO of a Schwarzschild black hole. Predictably, 
as the inclination increases, we see that both the abso- 
lute and relative amplitudes of the higher harmonics in- 
crease, almost to the limit of a periodic delta-function 
when / 90°. Again, we find very little dependence of 
the harmonic structure on hot spot size or shape, as long 
as the total emission of the spot relative to the disk is con- 
stant. Here we have normalized the rms amplitudes to the 
background flux from the disk with a hot spot overbright- 
ness of 100%. For a signal /(f) with Fourier components 

/(f) = a„cos{27tnt), (1) 

n=0 

we define the rms amplitude fl„(rms) in each mode n > 

as ^ 

fl„(rms) = — (2) 
V2flo 

With this normalization, the total rms can be conve- 
niently written 

rms = J'YaKims). (3) 



FITTING QPO DATA 

The high frequency QPOs observed with RXTE can be 
interpreted in the context of the hot spot model by match- 
ing the model parameters to the features of the QPO 



power spectrum, particularly the locations, amplitudes, 
and widths of the peaks. This problem of fitting a power 
spectrum to the parameters of a multi-dimensional model 
is analogous to the spectacular successes that have re- 
cently been made in observing and interpreting the per- 
turbations in the cosmological microwave background 
(CMB). Like the CMB power spectrum, the features of 
the QPO power spectrum often can be isolated and fit- 
ted to one or two of the model parameters at a time. We 
approach this problem in the systematic step-by-step pro- 
cedure outlined below, using the example of the Type A 
QPOs fromXTE J1550-564 |4] as case study. 

First, we match the frequencies of the QPO peaks 
to the coordinate frequencies of specific geodesic or- 
bits around a black hole in order to determine the black 
hole mass and angular momentum. The 2:3 frequency 
commensurability of the two major peaks (and possi- 
bly 2:3:4 in some observations) directs our attention to 
orbits with radial and azimuthal frequencies with a ra- 
tio of Vr : V)j, —1:3. The choice of the 1:3 ratio instead 
of 2:3 gives the appropriate relative power in the two 
peaks-a major peak at the fundamental azimuthal or- 
bital frequency with additional power in the sidebands 
at V = Vif, ±Vr- Conversely, a ratio of 2:3 in the coor- 
dinate frequencies would actually give a 1:3:5 ratio in 
the major power spectrum peaks. For any value of the 
black hole spin Q<a/M < 1 , there exists a special radius 
where nearly circular orbits will have the appropriate 1:3 
ratio in coordinate frequencies. Then by scaling the black 
hole mass, we can match the actual peak locations to the 
data, giving a one-dimensional degeneracy in mass-spin 
parameter space. 

This degeneracy could be broken by identifying the 
low frequency QPOs with the Lens-Thirring precessional 
frequency at the same special radius. However, this as- 
sociation is still quite speculative and has difficulty ex- 
plaining the large rms amplitude of the LFQPOs. For 
the present, we will constrain our model to a region of 
this degeneracy limited by the independent black hole 
mass determination of radial velocity measurements, giv- 
ing 8.5 <M/M0 < 11.5 1 13]. For the QPO peaks at 184 
Hz and 276 Hz in XTE J1550, we choose M = 1O.3M0 
and a/M — 0.5, giving a geodesic orbital radius of ro = 
4.89M, somewhat outside the ISCO at Risco = 4.23M. 

Next, we match the total rms amplitudes in the two 
major peaks by selecting an average hot spot overbright- 
ness of 100% and radius /?spot — 0.35M. The disk incli- 
nation is independently determined to be / « 70° 1 1 3] . 
As explained in the previous section, this relatively high 
inclination should produce significant power in many of 
the higher harmonics at integer multiples of the funda- 
mental frequency v = nv^j,. These harmonics lie well 
within the sensitivity range of RXTE, yet are not detected 
with such high amplitudes. Instead, we find that by shear- 
ing the hot spot into an arc of finite length in azimuth. 



the power at higher frequencies is strongly damped, due 
to the countering of the lighthouse effect when the emis- 
sion region is spread out over a larger portion of the orbit 
15]. An arc length of 180° gives the appropriate amount 
of power around 184 Hz and 276 Hz while also success- 
fully decreasing the power at higher frequencies. 

The relative power in the two major peaks is primarily 
determined by the eccentricity of the geodesic orbit, as 
well as the arc length of the hot spot emission region. For 
the Type A QPOs, even the modest eccentricity of 0.11 
is sufficient to produce significant power in the "side- 
band" at 184 Hz. Even the Type B QPOs, where the 
peak at 184 Hz is stronger than the fundamental peak at 
276 Hz, can be explained easily in the context of the hot 
spot model by shearing the arc almost to a ring-shaped 
perturbation with arc length 330°. For both cases, it is 
not obvious why there should be such a small amount of 
predicted power at + Vr =368 Hz, further advocating 
the necessity of a full ray-tracing algorithm, as opposed 
to simple analysis of the coordinate frequencies. 

Up to this point, all of the ingredients of the hot spot 
model combine to produce a power spectrum with a se- 
ries of discrete peaks describing a purely periodic light 
curve. In order to understand the significant broaden- 
ing of the peaks (and thus the term ^MfliZ-periodic), we 
generalize the hot spot model by replacing a single spot 
orbiting indefinitely on a single trajectory by a collec- 
tion of hot spots that are continually created and de- 
stroyed with random phase along their orbits. Each hot 
spot lifetime t has an exponential distribution function 
f{t) oc exp(— f/r) with a probability P = dt /T of being 
"destroyed" during each time step dt. When one hot spot 
is destroyed, a new hot spot is created with the same ro 
but with random orbital phase. The effect is to spread 
each delta-function feature in the periodic power spec- 
trum into a Lorentzian peak with FWHM °^ l/T. With 
this multiple hot spot model, the power spectrum from 
XTE J 1550 is closely matched by setting the typical life- 
time r = 15 ms, or approximately four orbits in azimuth. 

The summary of these best-fit results is shown in Ta- 
ble 1, and the corresponding power spectrum (calculated 
from a Monte-Carlo simulation of the light curve over 
many hot spot lifetimes) is shown in Figure|2 The solid 
curve is the best fit of a collection of Lorentzians to the 
simulated data, included a small background counting 
noise level with a Poisson distribution. Again, the anal- 
ogy to the CMB data is evident, and just like the with the 
CMB, each successive advance in instrumentation brings 
us closer to a precise determination of the model param- 
eters. However, for the QPOs, in addition to the power 
spectrum, there is also significant information in the light 
curve phase, if we can achieve the necessary sensitivity 
to resolve the signal in the time domain. This phase in- 
formation promises to ultimately help in distinguishing 
between the various proposed explanations of the cause 
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FIGURE 2. Simulated QPO power spectrum for a set of 
model parameters described in Table 1 , selected to fit the spec- 
tra of Type A QPOs from XTE J 1550-564. The solid line is a 
fit to the simulated data with a series of Lorentzian peaks. 

TABLE 1. Hot spot model param- 
eters for Type A HFQPOs from XTE 
J1550-564 



Parameter 


Value 


BH mass 


10.3Mp, 


BH spin a/M 


0.5 


disk inclination 


70° 


geodesic ro 


4.89M 


^spot 


0.35M 


eccentricity 


0.11 


arc length 


180° 


overbrightness 


100% 


hot spot Lifetime T 


15 ms 


(4 orbits) 




Lorentzian FWHM 


45 Hz 



of the black hole QPOs. 

We have also investigated the effect of various dis- 
tributions in orbital eccentricity, overbrightness, and the 
number of simultaneous hot spots at any time. Interest- 
ingly, we find that the QPO power spectrum is depen- 
dent only on the average eccentricity and overbrightness, 
and is not sensitive to the actual shape of the distribution. 
Similarly, the power spectrum is independent of the num- 
ber of different hot spots in existence at any given time, 
sensitive only to the average lifetime of each hot spot. 

One major remaining issue with the hot spot model is 
the preferred location of the geodesic that gives rise to 
1:3 coordinate frequencies. Why should the orbital fre- 
quencies favor integer ratios, and why should the pre- 
ferred ratio be 1:3 and not 1:2 or 1:4? It is possible that 
detailed radiation-hydrodynamic calculations with full 
general relativity will be required to answer this question. 
Perhaps the non-circular orbits can only survive along 
closed orbits such as these to somehow avoid destruc- 
tive intersections. For now, we are forced to leave this as 



an open question unanswered by the geodesic hot spot 
model. 

Yet with the computational framework developed 
above, these questions can be answered by modeling a 
whole collection of hot spots and arcs continually form- 
ing and evolving in shape and emissivity. The particular 
physical parameters for these hot spots can be derived 
from published MHD calculations such as |12|. With 
the basic ray-tracing and radiation transport methods in 
place, it should be possible to use our code as a "post- 
processor" to analyze the results of these 3D simulations 
to simulate X-ray light curves and spectra from a realistic 
accretion disk. 

A final piece of the black hole QPO puzzle is the spec- 
tral behavior of the source during outburst. Remillard et 
al. |4] have shown there is a large range of X-ray fluxes 
with different relative contributions from a power-law 
component and a disk bolometric component. These rel- 
ative and absolute fluxes seem to be correlated with the 
amount of power in both the LFQPOs and the HFQPOs. 
Future work on 3-dimensional disks and more detailed 
radiation transfer models should give us important in- 
sights into understanding this spectral behavior and its 
relation to the QPO power. 
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